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Abstract 

We present the currently most accurate prediction of Z pole observables such as 
sin^ ^cflf, Tzi Rb, Ri, and cr^^^^ in the Minimal Supersymmetric Standard Model (MSSM). 
We take into account the complete one-loop results including the full complex phase 
dependence, all available MSSM two-loop corrections as well as the full SM results. We 
furthermore include higher-order corrections in the MSSM Higgs boson sector, entering 
via virtual Higgs boson contributions. For r{Z XiXi) we present a full one-loop 
calculation. We analyse the impact of the different sectors of the MSSM with partic- 
ular emphasis on the effects of the complex phases. The predictions for the Z boson 
observables and Myy are compared with the current experimental values. Furthermore 
we provide an estimate of the remaining higher-order uncertainties in the prediction of 
sin^ ^efr- 
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1 Introduction 



Z boson physics is well established as a cornerstone of the Standard Model (SM) [1-3]. 
Many (pseudo-) observables [4] have been measured with high accuracy using the processes 
(mediated at lowest order by photon and Z boson exchange) 

e+e- ^ //, j + e, (1) 

at LEP and SLD with a center of mass energy -y/s Mz- In particular these arc the 
effective leptonic weak mixing angle at the Z boson resonance, sin^ ^cffj Z boson decay widths 
to SM fermions, V{Z — > //), the invisible width, Finv, the total width, F^, the ratios of 
partial widths, Ri and forward-backward and left-right asymmetries, A^-q and Alr, and 
the total hadronic cross section, o\^. Together with the measurement of the mass of the 
W boson, Mw, and the mass of the top quark, mj, the Z pole observables have been used to 
constrain indirectly the SM Higgs boson mass, M^'^, the last free parameter of the model, 
yielding Mf^ = 76l^| GeV with an upper hmit of Mfj^ < 144 GeV at the 95% C.L. [1, 
2]. The precision observables are also very powerful for testing models beyond the SM. In 
particular the Minimal Supersymmetric Standard Model (MSSM) [5] has been investigated, 
see Ref. [6] for a review. Performing fits in constrained SUSY models a certain preference 
for not too heavy SUSY particles has been found [7-12]. The prospective improvements in 
the experimental accuracies, in particular at the ILC with GigaZ option, will provide a high 
sensitivity to deviations both from the SM and the MSSM. 

In order to fully exploit the high-precision measurements, the theoretical uncertainty in 
the predictions of the (pseudo-) observables should be sufficiently smaller than the experimen- 
tal errors. Within the SM the complete one-loop and two-loop results [13-16] as well as lead- 
ing higher-order contributions [17-21] are available for sin^ ^eflf- For the leptonic and hadronic 
Z widths partial results for process-specific two- loop corrections are known [4, 22-24] . 

The theoretical evaluation of the Z pole observables within the MSSM is not as advanced 
as in the SM. So far, the one-loop contributions have been evaluated completely, restricted 
however to the special case of vanishing complex phases (contributions to the p parameter 
with non-vanishing complex phases in the scalar top and bottom mass matrices have been 
considered in Ref. [25]). At the two-loop level, the leading 0{aas) corrections [26] and the 
leading electroweak corrections of 0{al), O{ciit0ib), C'(q;^) to Ap have been obtained [27, 

28] {at and at are defined in terms of the Yukawa couplings yf as af — ^). Going beyond 
the minimal SUSY model and allowing for non-minimal fiavor violation the leading one-loop 
contributions to Ap are known [29]. 

In order to confront the predictions of supersymmetry (SUSY) with the electroweak pre- 
cision data and to derive constraints on the supersymmetric parameters, it is desirable to 
achieve the same level of accuracy for the SUSY predictions as for the SM. In this paper we 
present complete one-loop results for the Z boson observables in the MSSM with complex 
parameters, taking into account the full complex phase dependence. Besides the decays of 
the Z boson into quarks and leptons, we also provide a complete one-loop result for the 
partial width of the Z boson decay into the lightest neutralino, T{Z — > XiXi)^ which is the 
first one-loop result for this decay in the full MSSM. If the Z boson decay into the lightest 
neutralino is kinematically possible it contributes to the invisible width of the Z boson. We 
combine our new one-loop results with the full set of available higher-order corrections in the 
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MSSM. In order to recover the state-of-the-art SM results in the decouphng hmit where all 
supersymmetric particles are heavy, we consistently incorporate also those SM-type higher- 
order corrections which go beyond the results obtained within the MSSM so far. In this way 
we provide the currently most complete results for the Z boson observables in the MSSM 
(for the corresponding results for the W boson mass, M^-, see Ref. [30]). A public computer 
code based on our result for the electroweak precision observables (EWPO), i.e. the Z pole 
observables and Mw-, is in preparation [31]. It provides predictions for the EWPO in terms of 
the low-energy parameters of the MSSM, which can be freely chosen as independent inputs. 

We analyse the numerical results for the EWPO for various MSSM scenarios, such as 
SPS benchmark scenarios [32], scenarios with heavy scalar masses [33,34] and the CPX sce- 
nario [35]. The dependence of the results for the EWPO on the complex phases is investigated. 
For sin^ ^eff, showing the largest sensitivity of the Z pole observables to the SUSY loop cor- 
rections, we provide an estimate of the remaining theoretical uncertainties from unknown 
higher-order corrections. 

The rest of the paper is organised as follows: In Sect. 2 we introduce our notations and 
conventions. The Z pole observables are discussed in Sect. 3. Details about the evaluation of 
the higher-order corrections are given in Sect. 4. In Sect. 5 we present our numerical analysis. 
The estimate of the theory uncertainties in the prediction for sin^ ^eff from unknown higher- 
order corrections is given in Sect. 6. We conclude with Sect. 7. 

2 Notations and conventions 

In the complex MSSM, the Z pole observables depend on all the free parameters of the model, 
such as SUSY particle masses, mixing angles and couplings. In order to be self-contained we 
list in this section the notations and conventions used in our calculations. We briefly describe 
the relevant quantities in the sfermion, the chargino/neutralino, and the Higgs boson sector 
of the MSSM. 

2.1 Sfermions 

The mass matrix for the two sfermions of a given flavour, in the /l, basis, is given by 



where {cot /3 , tan/?} applies for up- and down- type sfermions, respectively, and tan/3 is 
the ratio of the two vacuum expectation values of the two Higgs doublets (see Sect. 2.3), 
tan/3 = V2IV1. We have furthermore used = sin^ = 1 — M^/M|. In the Higgs and 
scalar fermion sector of the complex MSSM, Nf+1 phases are present, one for each Af and 




(2) 



with 



Ml = Ml + Ml cos 2/3 (// - g/4) , 
M| = M|, + M| cos 2/3 g/4, 
Xf ^ Af - 11* {col P , tan /3}, 



(3) 
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one for i.e. Nf + 1 new parameters appear in comparison to the MSSM with only real 
parameters. As an abbreviation, 

(f>Af = arg(A^) (4) 

will be used. As an independent parameter one can trade (pAf for (f)Xf = arg(X/). The 
sfermion mass eigenstates are obtained by the transformation 



with a unitary matrix Uj. The mass eigenvalues are given by 

1 
2 



/1, 2 



m 



f 



Ml + Ml^ J {Ml - MIY + Am)\Xf\ 



and are independent of the phase oi Xj. 



(5) 



(6) 



2.2 Gaugino sector 

The physical masses of the charginos are determined by the matrix 

M2 \^s\n(3Mw 
cos /3 Mw A* 



X 



(7) 



which contains the soft breaking term M-i and the Higgsino mass term /x, both of which may 
have complex values in the complex MSSM. Their complex phases are denoted by 



!)M2 = arg (M2) and 0^ = arg (/i) . 



(8) 



The physical masses are denoted as "^-^^ and are obtained by applying the diagonalisation 



matrices U^± and Vj^± 



(9) 



The situation is similar for the neutralino masses, which can be calculated from the mass 
matrix (s^ = sin c^ — cos 6*^) 



/Ml -MzSwCos/3 MzSwSin/3 ^^ 

M2 M^Cw cos/? -M^Cw sin/? 

—M^s„ cos/? M^c„ cos/? —\i 

\ M^Sw sin/? —M^Cw sin/? — / 



(10) 



This symmetric matrix contains the additional complex soft-breaking parameter Mi, where 
the complex phase of Mi is given by 



0Mi = arg (Ml) 



(11) 



The physical masses are denoted as m 



dure using the matrix N^o 



Xl,2,3,4 



and are obtained in a diagonalisation proce- 



N^oYNto = diag (|m^o,m^o,m^o,m^o) . 



;i2) 
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The gluino enters the predictions for the hadronic decays of the Z boson (and accordingly 
the total Z width) at the one-loop level, while for the EWPO with leptons in the final state 
it enters only at 0{aa^. The soft-breaking gluino mass parameter M3 is in general complex, 

M3 = |M3|e^'^^ (13) 

and the gluino mass is given by mg = IM3 1 . The phase can be absorbed by a redefinition of the 
gluino Majorana spinor such that it appears only in the gluino couplings but not in the mass 
term. In our calculation of the EWPO below we will incorporate the full phase dependence 
of the complex parameters at the one-loop level, while we neglect the explicit dependence 
on the complex phases beyond the one-loop order. Accordingly, we incorporate the gluino 
phase appearing in eq. (13) into our predictions for the hadronic Z observables (and also in 
Higgs-sector corrections associated with the bottom Yukawa coupling, see below), while we 
treat the two-loop corrections to the EWPO in the approximation of vanishing gluino phase. 



2.3 Higgs bosons 



CT'-violating phases can have an important impact on the Higgs sector of the MSSM with 
complex parameters [36-40]. In our (one- loop) analysis of the complex phase dependence 
of the predictions for the Z pole observables we take into account also the potentially large 
effects of CP-mixing in the complex MSSM Higgs sector [40,41], although formally CP- 
violating contributions in the complex MSSM Higgs sector enter the Z pole observables only 
at the two-loop level. We consistently incorporate the higher-order corrected Higgs states in 
our analytical results for the Z pole observables and the W boson mass. 

Once higher-order terms are included in the complex MSSM [36], Higgs fields are mixtures 
of the CP-even, h and H, and the CP-odd states, A and G. The Higgs boson propagator 
matrix receives contributions from the Higgs boson self energies. In the propagator matrix 
the mixing between the Higgs fields h,H,A and the Goldstone boson G is of sub-leading 
two- loop order and can therefore safely be neglected [40,42]. The three mass eigenvalues of 
the remaining (3 x 3) matrix, 

Mh, < < Mh„ (14) 

corresponding to the mass eigenstates hi, h2, h^, are then determined by the propagator 
poles [40]. 

It is sometimes convenient for phenomenological analyses to introduce effective couplings 
that incorporate leading higher-order effects. There is no unique procedure how to define 
these effective couplings, and care has to be taken not to spoil gauge or unitarity cancellations 
by a partial inclusion of higher-order contributions. One possibility is to use the "p^ = 0" 
approximation (see Ref. [40] for details), i.e. to evaluate all Higgs boson self energies at zero 
external momentum, which is equivalent to the effective potential approach. In this way a 
unitary effective mixing matrix Ucff can be defined (ensuring decoupling to the SM for heavy 
SUSY particles), transforming in this approximation the lowest-order states h, H, A into the 
mass eigenstates hi, h2, h^. 



( hi\ 



eff 



H 



(15) 
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The elements of the effective mixing matrix Ues can be interpreted as effective couplings 
of the Higgs bosons, incorporating leading higher-order corrections from Higgs boson self 
energies into the Higgs boson couplings [40]. In our numerical calculation the Higgs sector 
parameters are evaluated with the help of the program FeynHiggs [40,43,44]. 

Another numerically important correction appears in the relation between the bottom- 
quark mass and the bottom Yukawa coupling, yi,. The leading tan /3-enhanced contributions 
to the relation arise from one-loop contributions with gluino-sbottom and chargino stop 
loops. We include the leading effects via the quantity Aj, [45] (see also Refs. [46-48]). The 
Afe corrections can affect for instance the prediction for T{Z —>■ bb), where yf, enters at the 
one-loop level. Thus the A{, corrections are a two-loop effect, which, however, can in principle 
have a noticable impact. 

Numerically the correction expressed by to the relation between the bottom-quark 
mass and the bottom Yukawa coupling is usually by far the dominant part of the contributions 
from the sbottom sector (see also Refs. [49, 50]). In the limit of Mp, Mp, » rrit and tan /3 » 1, 
Afe is given by [45] 

= Y^M^^f tan/3 x I{mi^,mi^,m-g) + ^ i^* tan/3 x I{mi^,mi^, |/x|) . (16) 
The function / is defined as 

For the numerical evaluation we use the implementation of A^ and the corresponding Higgs 
couplings in FeynHiggs. 



3 Description of the Z boson resonance 
3.1 Effective coupling approach 

e^e~ collisions at the Z boson resonance are commonly described in an effective coupling 
approach. These effective couplings are subsequently used to define the so-called Z resonance 
pseudo observables. In the following we will mostly use the conventions given in Ref. [4]. At 
Born level the matrix element of the process in eq. (1), depicted in Fig. 1, is given by 

A^Born OC -{QeQf {la ® l") + X Wv ,{0)91 (7a ® 7^) " ^(0)4(0) (7a ® 7^75) 



(18) 



-5a,(0)^((0) (7a75 ® 7") + ^a,(0)4(0) (7^75 ® 

with the propagator x defined as^ 

^ " s-M^ + isTz/Mz' ^^^^ 

"'^Thc form of the propagator in cq. (19) corresponds to a Breit-Wigner function with a running width, 
which is the convention normally adopted in the experimental determination of the gauge-boson masses. 
It should be noted, however, that from two-loop order on the Z boson propagator has a complex pole. 
Expanding around the complex pole and defining the mass according to the (gauge-invariant) real part of 
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In eq. (19) s denotes the center-of-mass energy, Mz is the Z boson mass and Fz its width. 
We use the shorthand notation 



{Aa ® 5") = [VeA^Ue] X [UfB"Vf] . (20) 

Furthermore the lowest-order vector and axial vector couplings of a fermion / 

f _ 4 - 2Q/4 

4(0) - ^> (21) 

were introduced in eq. (18). Qf denotes the charge of the fermion / as fraction of the 
elementary charge e, // is its third weak isospin component. The simple structure of the 
matrix element in eq. (18) easily allows the identification of QED contributions due to photon 
exchange graphs (first term) and electroweak Z boson exchange graphs (terms in square 
brackets) . 




The inclusion of higher-order corrections in general leads to a modification of this simple 
structure. At the one-loop level one can easily classify the corrections into QED, QCD 
and electroweak corrections. Supplementing the Born-level diagrams (see Fig. 1) with an 
additional photon line gives the one-loop QED corrections. They form a gauge-invariant, UV- 
finite, but IR-divergent subset of loop corrections to the process e+e~ — > //. IR finiteness is 
obtained by taking real Bremsstrahlung contributions into account {e'^e~ — > ff + j). This 
implies that the QED corrections are dependent on the experimental setup. The situation 
is similar for the QCD corrections. At one-loop order these corrections occur for quark pair 
production (e^e~ — > qq) and contain corrections due to virtual gluon and gluino exchange, 
as well as real gluon emission, the latter again leading to a dependence on the experimental 
setup. Beyond one-loop order a clean distinction between QCD and electroweak corrections 
is no longer possible. Mixed 0{aas) corrections appear for instance as gluon- and gluino- 
exchange contributions in virtual quark loops of the W and Z boson propagators. In the 
following we will explicitly indicate final-state QCD corrections and initial- and final-state 

the complex pole leads to a Breit-Wigner paramctrisation of the resonance line shape with a constant decay 
width (see Ref. [51] for details). The (numerically sizable) difference between the two mass definitions needs 
to be properly taken into account when incorporating higher-order corrections obtained with the fixed-width 
parametrisation. 
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QED corrections where relevant, while we use a common notation for all other contributions. 
The self energy, vertex and box contributions can be expressed in terms of complex form 
factors T^j . These form factors depend on the Mandelstam variables s and t, where the 
t dependence enters only via box contributions. 

At the Z boson resonance those contributions that are not enhanced by a resonant Z prop- 
agator are relatively small. In particular, box diagrams contribute only a fraction of less than 
10~^ at the one-loop level. It is therefore convenient for describing physics at the Z boson 
resonance to treat the non-resonant higher-order corrections separately as part of a "de- 
convolution" procedure (see the discussion in Sect. 3.2 below) and to express the dominant 
contributions in terms of a Born- type matrix element. The effective couplings in this matrix 
element are given by the form factors T^- = V, A) in the approximation where higher- 
order non-resonant contributions are neglected, so that they only depend on the Mandelstam 
variable s, 



M^''^-^ a{s) (7a ® 7") + X ^vvis) {la ® 7") - ^v'a{s) {la ® 7^75) 



1 r 



s y 



-^%{s) (7a75 ® 7") + ^li(^) (7a75 ® 7^75) ] }• (22) 



As the only relevant contribution to the form factor of the photon-exchange part, the running 
QED coupling a{s) is kept. 

The factorization is the result of a variety of approximations that is valid at the Z reso- 
nance to the accuracy needed (see Ref. [4] for a more detailed discussion). At this level, the 
form factors T^- in eq. (22) factorise into contributions from the production and the decay. 
The final step in the "Z pole approximation" is to set s — M| in the effective matrix element 
of eq. (22). This yields for the form factors 

Ttl{s = Ml) = glsl (23) 
where the couplings g^^y^-^ have the loop expansion 

1 + 4(1) +4(2) + ■■■ 



f — f 

9v — 9v,{o) 

f — f 

9 a — 9a,{0) 



1 + 4(1) + 4(2) + •••!• (24) 



We use indices {V, ^4} in capitals to label generic vector and axial vector couplings which 
contain all higher order terms. In contrast, lower case {v, a} is used for couplings of a specific 
loop order. The couplings g'^^ ^^^^ in the above equation are thus the lowest-order couplings 

from eq. (21). The terms 4'ua}{(i)(2) } represent loop corrections of order {(1), (2), . . . }. 

Note that, working in our conventions, the Born level couplings 9^^^^^} factored out in 

eq. (24) and are hence not contained in the higher order terms g'^^^y }• 

The Z boson decay, which almost entirely proceeds through the decay into two fermions 
(the Born level graph is depicted in Fig. 2) 

Z ^ //, (25) 
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is with these prerequisites of the form 

^^e«, (26) 

where Uf,Vf are the Dirac spinors of the fermion anti-fermion pair and is the polarisation 
vector of the Z boson. 



9v - 9^5 




Figure 2: Z boson decay to fermions at Born leveL 



3.2 Observables at the Z boson resonance 

The Z pole (pseudo-) observables determined from the measurements at the Z boson reso- 
nance in e'^e~ annihilation allow to perform precision tests of the SM [4] and the MSSM [6]. 
The most prominent Z pole observables, some of which are interrelated (see the discussion 
in the subsections below), can be classified as follows: 

(a) inclusive quantities: 

— the partial leptonic and hadronic decay widths F/, 

— the total decay width F^, 

— the hadronic peak cross section cr^^^, 

— the ratio of the hadronic to the electronic decay width of the Z boson, Ri, 

— the ratio of the partial decay width for Z — > cc {bb) to the hadronic width, i?c(6)- 

(b) asymmetries and effective fermionic weak mixing angles: 

f 

— the forward-backward pole asymmetries A-^^, 

Of 

— the left-right pole asymmetries ^l'^^, 

— the effective fermionic weak mixing angles sin^^^g. 

The quantities that can be measured directly in collider experiments are energy dependent 
cross sections (T-^{s), left-right A(^^{s) and forward-backward asymmetries Apg(s) of the pro- 
cesses e+e~ — >• (7, Z) — > ff{n'j) / ff{ng) and the corresponding Bhabha reactions for f = e. 
Applying deconvolution (unfolding) procedures to these quantities, pseudo-observables at the 
Z resonance such as masses, partial widths, pole asymmetries and effective mixing angles can 
be obtained. The deconvolution procedure is performed under certain assumptions, in the 
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present case the Z pole approximation, with the support of sophisticated programs such as 
ZFitter [52] or TopaZO [53]. This means in particular that by applying the Z pole approx- 
imations none of the information from the theory side is lost, but is already accounted for 
when the pseudo-observables are calculated from the raw data. 



3.2.1 Effective mixing angles 



The effective electroweak mixing angles of fermions /, in particular leptons, are of great 
relevance for testing the SM and its extensions. They have been measured with high accu- 
racy, while the theory predictions sensitively depend on the model under consideration. The 
effective electroweak mixing angle of a fermion / is defined as 



4|Q/|sin'^fg := 1-Re 



(27) 



This definition allows to absorb all higher-order corrections into a quantity Ak/ (similar to 
the case of the quantity Ar in the context of the muon decay, cf. Ref. [30]), 



sin' 




1 - Re 



9± 



4Re [k/] = 4(1 + Re [Ak/]), 



(28) 



which expanded up to two-loop order reads 

^«^/ = (^((1) - <(!)) + (4(2) - 4(2) + - <(!))) 



+ 



(29) 



The theoretical prediction for the quantity Ak/ can be obtained in different models, depend- 
ing on all free parameters of the model through virtual corrections. 



Ak/ = Afi;/(MH/, Mz, mt, a, a^, . . . , X), 



(30) 



where 



X 



Mf^ (SM), 



X = Mh, Mh, Ma, Mh±, tanp, Mf, Af,m^o,±, 
Furthermore the prefactor 



(MSSM). 



1 - 



Mi 



w 



1 - 



M^(Mz,mt,a,as,...,X) 



Ml 



(31) 



itself depends on the experimental input M^, mt, a, c^s, • • • and the unknown model param- 
eters X. This is due to the fact that the theory prediction for Mw is employed, which is 
obtained from muon decay by trading the precisely measured Fermi constant as input for 
Mw (see Ref. [30] for the most up-to-date prediction for My^ in the MSSM) . 
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3.2.2 Asymmetries and asymmetry parameters 



Though directly related to the effective mixing angles, there are further observables which 
are often referred to in the literature. For completeness we include them in our discussion. 

We start with the asymmetry parameters Af. They are by definition given in terms of 
the real parts of the ratios of the effective vector and axial vector couplings, in contrast to 
the partial widths that include the full complex effective couplings. 

Re 4 

Af = 2 (32) 

l+[Re4]' 

The asymmetry parameters Aj arc obviously directly correlated with the effective mixing 
angles defined in eq. (27). They can be used to calculate the forward-backward pole asym- 
metries Ap'l and the left- right pole asymmetries A^^, 

A'^i = ^^AeAf, (33) 
A^i = Ae. (34) 

Note that the Ap^ are theoretical, to a certain level artificial, quantities (pseudo-observables 
in the terminology of Ref. [4]). Besides the unfolding from the QED corrections, also small 
effects from j-Z interference and pure 7 exchange have to be taken into account, as men- 
tioned in the beginning of Sect. 3.2. Their use is justified by their dominance at the Z peak, 
by their simple parametrizations in terms of the effective couplings resp. mixing angles, and 
by the practically model-independent correction terms relating them to realistic observables. 
For a detailed discussion, see also Ref. 



3.2.3 Decay widths 

The partial decay width for the decay into fermions / is given by (cf. Ref. [4]) 



/ 
9v 



R 



f 
9 a 



Ri 



(35) 



where the radiation factors Ryj^ describe the final state QED and QCD interactions and 
account for the fermion masses ruf, see eq. (65) below. The latter contributions are of 
particular relevance for the decay into bottom-quarks [54]. 

Equivalently the partial decay width can also be expressed as (cf. Ref. [55]) 



= TV/f \pf\ (4(// - 2Qfsl \Kf\fRi + R{) , 



with 



(36) 



(37) 



24V27r' 

The factor normalises the overall matrix element to Z boson decay. The quantity Kf is 
a measure for the relative coupling strength of the effective vector and axial vector couplings. 
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see Sect. 3.2.1. The effective couphngs gf^^ used in eq. (35) and the quantities p/ and kj 
from eq. (36) are related via 

l-MQf\sl^f=4- (38) 
9a 

It is convenient to express pf and in terms of universal corrections (that are independent 
of the fermion species) and a non- universal part (depending on the fermion species), 

P/ = 1 + Apf = 1 + Apuniv + Ap/,non-univ, 

K/ = 1 + Ak/ = 1 + A/«univ + Are/^non-univ (39) 

The universal part arises from vector-boson self energies and fermion-independent countert- 
erms, and the non-universal part from vertex corrections and fermion-dependent counterterms 
(see e.g. eq. (56)). The leading contribution to the universal corrections arises from mass 
splitting in isospin doublets, like an isospin doublet of a given family of SM fermions or 
MSSM sfermions, entering via the quantity Ap [56] according to 

Apuniv = Ap . . . , 

A«univ = %Ap+... , (40) 

where 

contains the transverse parts of the Z and W boson self energies evaluated at vanishing 
momentum-squared. In the approximation where all fermion masses except the top-quark 
mass are neglected (m^ ^ rrib, rrir, ■ ■ ■), Ap in the SM, at the on- loop level, is given by the 
simple expression 

2 

Expressing the weak coupling in terms of the Fermi constant rather than the fine-structure 
constant a yields 

where the bar is used to distinguish the two parametrisations. Leading reducible two-loop 
terms can be obtained by performing the substitutions [57] 

2 2 

«:^ = l + %Ap ,i + ^Ap + ... . (44) 
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It should be noted that the substitutions in eq. (44) replace Ap (parametrised in terms of 
a) by Ap (parametrised in terms of G^. This introduces additional terms at higher orders 
according to 

Ap(a) ^ Ap(G'^) = (1 + Ar)Ap(a), (45) 

which follows from the explicit expressions given in eq. (42) and cq. (43). 

It was already mentioned that the Z predominantly decays into fcrmion pairs. Contribu- 
tions from other decay modes such as Z — >■ ^^'f,ggg,bbbb are insignificant [58]. In the SM 
the total width Fz is thus given as the sum of the leptonic width Fi and the hadronic width 

rl^^r^ + r^, (46) 

where 

r/ = rinv + re + r^ + r„ (47) 

and 

= r„ + + Te + r, + f^ (48) 

The invisible width Finv in the SM used in eq. (47) arises from the partial widths of 
Z boson decays into neutrinos, 

ri„v = r,, + r,^ + r,,. (49) 

In the MSSM further contributions to the invisible width can arise from the decay into the 
hghtest neutralino (if this particle does not decay further within the detector) . We therefore 
include the process Z — > XiXi into our computations, where x° is the lightest neutralino. 
This decay channel of the Z boson is open if m^o < Mz/2. The total Z boson width in the 
MSSM thus reads 

rz = r; + r;, + r^o, (50) 

where the partial width F^o is in analogy to eqs. (35), (36) given by 



9 a 



'nf, r^o = ro|p^o|i?f, (51) 



x" x" x^ 

with the radiation factor R/ and the effective axial vector coupling specified below, g^^ is 

related to p^o via eq. (38), with f ^ xl- As neutralinos are described in terms of Majorana 

fermions, the coupling of Xi to the Z boson does not contain a vector part. 

3.2.4 Peak cross-sections and ratios of partial widths 

The cross-section for the process eq. (1) is given by 



-? = 12-^, (52) 



at the Z boson resonance, i.e. s = Mz- The hadronic peak cross section is obtained in 
complete analogy as 

- ''^Wl- '''' 
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the only difference being tliat all decay modes into hadronic final-state particles are consid- 
ered. Furthermore certain ratios of partial widths are often analysed, 



Ri = ^, (54) 

'- e 

Rb,c = ^, (55) 

i h 

with the partial leptonic (electronic) Fg and hadronic F^, Fc widths from eq. (47) and eq. (48), 
as well as the partial decay widths to charm Fc and bottom-quarks F^. 



4 Calculation of effective couplings and pseudo observ- 
ables 

4.1 Complete one-loop result in the complex MSSM 

In order to calculate the Z pole observables, one of the main tasks is to evaluate the effective 
couphngs g'^^Y A}- ^^^^ requires the computation of the Zff vertex graphs (see Figs. 17, 18, 
19 given in the Appendix) and the self energy, as well as the evaluation of the respective 
countcrterms. Expressed in terms of unrenormalised vertex graphs and on-shell counterterms 
the effective couplings expanded up to one-loop order are found to be 



si 



(a) 



s'v 



(a) 



9a,{Q) 
9a,{0) 



(1) 



^ ^ (vertex) f ^8e^5sl si - 4 



9a,{0) 



2cl 



2^^ 



(0) 



9v,{o) 



9vM 



1 + 4 



(1) 



(vertex); 
/ 

9v,{Q) 



6 S-tTiT 



2QfS 



+-SZ^^ + 

Qf (T?^{M^ 



2cl {li - 2Qfsl) 



iQf 1 



USZ{ + 6Z 



2ff,(o) 



9v,{o) 



(56) 



where (vertex)^^^ stands for the axial/vector part of the unrenormalised vertex graphs, 
schematically depicted in Figs. 17 and 18. The exphcit form of the on-shell counterterms 
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in eq. (56), expressed in terms of self-energies, can be found in Ref. [59] (we use the conven- 
tions given therein). They have the same structure as in the SM, but the self-energies have 
to augmented by the non-standard contributions, accordingly. The charge renormalisation 
counterterm, 5e, receives contributions of large logarithms from light fermions, giving rise 
to a numerically sizable shift in the fine-structure constant. Act = Attj;^^ + Aaicpt- The 
correction Ap, see eq. (41), enters via the counterterm to the weak mixing angle, (5s^. 

Computing the occurring Feynman graphs we include the complete set of MSSM one- 
loop diagrams, keeping the full complex phase dependence. All relevant Feynman graphs 
are calculated making use of the packages FeynArts [60] and FormCalc [61]. As regulari- 
sation scheme dimensional reduction [62] is used, which allows a mathematically consistent 
treatment of UV divergences in supersymmetric theories at the one-loop level. Our results 
are presented in the same conventions as in Ref. [30]. The higher-order corrections in the 
Higgs sector were implemented into our predictions for the Z pole observables and M^/ as 
described in Sect. 2.3. Concerning the light quarks and leptons, in our calculations we take 
into account the masses of the bottom-quark and the tau-lepton, while the masses of the 
other leptons and light quarks are neglected (except for the contributions giving rise to the 
shift Act in the fine structure constant, as discussed above). As a consequence of keeping 
a non-zero bottom-quark and tau-lepton mass, in addition to the terms given in cq. (22) 
the matrix elements in our calculations contain also contributions oc UT,biki)urfi{k2){ki^2^z) 
and oc Urfi{ki)'y^Urfi{k2) {ki^2^z) ■ We checked by explicit computation that these terms are 
numerically negligible. 

As mentioned above, the decay channel Z — XiXi can provide additional contributions 
to the invisible width of the Z boson and thus has to be taken into account. In order to 
obtain results at the same level of precision as for the leptonic and hadronic observables, we 
have performed for the first time a full MSSM one-loop calculation for the decay Z XiXi- 
At Born level the axial vector couplings are given by (there are no further contributions from 
vector couplings owing to the Majorana nature of the neutralino) 




(57) 




+4c^.^(vertex)f? + §^^^(A^i,3A^r,3 - N^,,Nl,) 




(58) 
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Here Nij{i,j = 1 . . .4) represent the in general complex entries of the neutralino diagoual- 
isation matrix N^o, see eq. (12). We use the notation (vertex)^" to denote the axial part 
of the vertex graphs displayed in Fig. 19. The counterterms of charge, electroweak mixing 
angle, and Z field renormalisation are again evaluated in the on-shell scheme. For the neu- 

tralino field renormalisation constants SZj^'' we employ the on-shell renormalisation detailed 
in Refs. [63,64]. The parameter tan/? is renormalised by imposing a DR condition as speci- 
fied in Refs. [65,66]. Analytic expressions for the neutralino field renormalisation constants 
in eq. (58) can be found in Ref . [63] . 



4.2 Incorporation of higher-order contributions 

Having computed the results for the effective couplings at the one-loop level (for the first time 
under consideration of the full complex parameter dependence, CV mixing in the Higgs sector, 
and resummed tan [3 enhanced Yukawa couplings) we now include the available higher-order 
contributions in the SM and the MSSM. As a result, we obtain the currently most accurate 
prediction of the Z pole observables in the MSSM. 



4.2.1 Combining SM and MSSM contributions 

As mentioned before, the theoretical evaluation of the Z pole observables in the SM is signif- 
icantly more advanced than in the MSSM. In order to obtain the most accurate predictions 
within the MSSM it is therefore useful to take all known SM corrections into account. This 
can be done by writing the MSSM prediction for a quantity x — gy^, pf, Kf, . . . as 

„MSSM _ ^SMI , ^MSSM-SM ^ ^SM I , ™SUSY /f-Q\ 

where x^^ is the prediction in the SM with the SM Higgs boson mass set to the hghtest 
MSSM Higgs boson mass, Mh^, and ^^^^''^"^^ = x^^^^ denotes the difference between the 
MSSM and the SM prediction. 

In order to obtain x^''^^^' according to eq. (59) we evaluate a;'*^^^'^"^'^ at the level of 
precision of the known MSSM corrections, while for x^^ we use the currently most advanced 
result in the SM including all known higher-order corrections. As a consequence, x^^ takes 
into account higher-order contributions which are only known for SM particles in the loop, but 
not for their superpartners (e.g. two-loop electroweak corrections to Ak beyond the leading 
Yukawa contributions). 

It is obvious that the incorporation of all known SM contributions according to eq. (59) 
is advantageous in the decoupling limit, where all superpartners are heavy and the Higgs 
sector becomes SM-likc. In this case the second term in eq. (59) goes to zero, so that the 
MSSM result approaches the SM result with Mf^ = Mh^. For lower values of the scale 
of supersymmetry the contribution from supersymmetric particles in the loop can be of 
comparable size as the known SM corrections. In view of the experimental bounds on the 
masses of the supersymmetric particles (and the fact that supersymmetry has to be broken), 
however, a complete cancellation between the SM and supersymmetric contributions is not 
expected. Furthermore, the leading Yukawa enhanced corrections with MSSM Higgs boson 
exchange were found to be very well approximated by the corresponding Yukawa terms in 
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the SM [27,28]. Therefore it seems appropriate to apply eq. (59) also for rather light SUSY 
particle spectra. 

4.2.2 Universal SM contributions beyond one-loop 

It is convenient to parametrise higher-order SM corrections to the Z pole observables in terms 
of the quantities and kj, defined in eq. (38). For the universal contribution Ar, entering 
via Miy, the complete two-loop result is available in the SM [51,67-69]. Further reducible 
SM higher-order contributions, appearing for instance as products of one-loop contributions 
in two- loop counterterms, can be incorporated with the help of substitutions like eq. (45). 
Beyond the two-loop order, irreducible higher-order corrections have been obtained for Ap, 
the leading universal contribution from the mass splitting in an isospin doublet (for two-loop 
contributions to Ap, see Refs. [70-73]). Higher-order QCD contributions to Ap in the SM arc 
known up to Oipto?^ [17]. The 0{ac?^ corrections to Ar are also known [18]. Electroweak 
three-loop contributions to Ap of C(mjG^) and mixed electroweak and QCD corrections of 
0{mfGf^as) have been obtained in Refs. [19,20]. Most recently even the full class of four- 
loop 0{aal) contributions to Ap became available [21] and has been included, although the 
corrections turned out to be rather small numerically. 

4.2.3 Universal MSSM two-loop contributions 

Within the MSSM, leading irreducible two-loop contributions to the Z pole observables arc 
only available as universal corrections to Ap in the approximation where complex phases are 
neglected. Reducible higher-order terms can be obtained in the same manner as in the SM. 

Leading irreducible SUSY QCD corrections of 0{aas) entering via the quantity Ap arise 
from the diagrams shown in Fig. 20. They involve both gluon and gluino exchange in (s)top- 
(s)bottom loops and were first evaluated in Ref. [26]. Besides the 0{aas) contributions, also 
the leading electroweak two-loop corrections of 0{a'^), 0{al) and 0{atab) to Ap have become 
available [28]. These two-loop Yukawa coupling contributions are due to MSSM Higgs and 
Higgsino exchange in (s)top-(s)bottom-loops, see Fig. 21. In Ref. [28] the dependence of the 
0{a^j^) corrections on the lightest MSSM Higgs boson mass, Af/^, was analysed. Formally, at 
this order the approximation Mh = would have to be employed. However, it was shown 
in Ref. [28] how a non-vanishing MSSM Higgs boson mass can be consistently taken into 
account, including higher-order corrections. Correspondingly we use the result of Ref. [28] 
for arbitrary and employ the code FeynHiggs [40,43,44] for the evaluation of the MSSM 
Higgs sector parameters. 

The final step is the inclusion of the complex MSSM parameters into the two-loop results. 
So far all generic two-loop results have been obtained for real input parameters. Following 
Ref. [30], we approximate the two-loop result of an observable Oz — Fj, F^, sin^ O^g, ... for 
a certain value of phase by a simple interpolation, based on the full phase dependence at 
the one-loop level and the known two- loop results for real parameters, O^^^(O), 0^^^(7r), 

0^(0) = on<l>) + [ono) - o^o)] X i^l^ 

+ [Ol'V) - Of{.)] X '-^^ . (60) 
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Here 0^{(j)) denotes the one- loop result, for which the full phase dependence is known. The 
factors involving cos0 ensure a smooth interpolation such that the known results O^^^(O), 
0|'"(7r) are recovered for vanishing complex phase. As a check the formula has been applied 
to the one-loop case. The numerical difference between the approximated and the full one- 
loop result was at the level of a typical SUSY two-loop contribution, which is expected for a 
pure one-loop result. 

4.3 Leptonic observables 

We now combine the various contributions discussed in the previous sections in order to 
obtain results for the leptonic observables at the Z boson resonance. 

4.3.1 Effective leptonic weak mixing angle sin^ 0eff 

In the SM the evaluation of higher-order corrections to Ani and sin^ ^eff is far advanced 
(here and in the following we use the notation sin^^gfr = sin^^gg). Recently the complete 
(9(q;^) contributions have become available. They include fermionic contributions [13,14], i.e. 
graphs with at least one closed fermion loop, as well as bosonic corrections [15, 16] without 
closed fermion loops. The 0{aas) and 0{aal) results can be found in Ref. [18] in terms 
of a leading top mass expansion. Universal three- and four-loop order corrections entering 
via Ap arc incorporated following the prescription in cq. (40). In total, the statc-of-thc-art 
expression for A/t^^^ can be decomposed into the following contributions (assuming lepton 
universality, we drop the index / in the following) 

A«SM ^ (A«;("))S^^ + (A«("'))S^^ + (A«(°"^))S^ 

+ (A«;("°^'))S^ + ^ (^Ap^"'*^>^^ + /^p^^^V) + Ap(°^"')) . (61) 

Accordingly, the SM prediction for sin^ ^eflf is given by 

Sin^ ^fff'laW+aas+aai+aaH'n-tG'jias+m^Gji = ("^w"^)^ (1 + Ak^^) . (62) 

As already mentioned in Sect. 3.2.1, it is crucial to use the theoretical prediction for My/ to 
calculate and Ak, employing in this way the more precisely measured experimental input 
parameter rather than . A simple parametrisation formula for (Ak*^" )^sm^ which 
approximates the full 0{a^) result to a precision even below the anticipated ILC/GigaZ 
errors, is given in Ref. [16]. In our analysis we use this result for the SM prediction of Ak as 
well as the currently most precise SM prediction for M\y [69]. 
Using eqs. (59) and (61), our result for Ak in the MSSM reads 

A« = A«:^ssM ^ A«:SM 1^^^^^^^ + (A«:H)Susy + ^ (^(Ap("«^))«usY ^ (^^(a^^susY^ _ (53) 

In this way the complete 0{a'^) contribution in the SM is incorporated, as well as the 
higher-order corrections which are either only known in the SM or solely for real MSSM 
parameters. Accordingly, our prediction for sin^ ^efr in th^ MSSM is given by 

Sin^ ^eflf |a+a2+aas+aai+aa3+^4(32c^^+^6G3 = S^(l Ak), (64) 
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where the most accurate MSSM prediction for Mw [30] is used in the calculation of 
and Ak. 

4.3.2 Leptonic decay widths 

In our calculations F;, the leptonic decay widths, are obtained in terms of the form factors 
Pi, Ki as described in Sect. 3.2.3. The form factor can be directly related to the leptonic 
mixing angles from cq. (64). The overall normalisation of the decay width, p^, is calculated 
from the effective one-loop couplings g\y\^ as described in Sect. 3.2.3. It is supplemented 
with the leading universal SM and MSSM corrections from Sect. 4.2 by applying eq. (44). 
Subleading corrections of 0{aas) are also available in the hterature [74]. The charge renor- 
malisation counterterm, as mentioned above, contains the numerically sizable shift in the 
fine-structure constant, Aa = A.a\^^^ + Acticpt. As these contributions are accounted for as 
part of 5'{y^}, they are consequently also incorporated into pi. In a final step the above 
results are supplemented with the radiation factors For leptonic final states these do 

not contain QCD corrections and are of the simple form (see, for example, Ref. [55]) 



with the running electromagnetic coupling constant a at the scale Mz- The inclusion of 
lepton masses mi yields numerically negligible effects. 

4.4 Hadronic observables 

The calculation of the hadronic observables proceeds in a manner very similar to the leptonic 
case. The discussion is therefore kept very short, addressing only some additional features in 
the hadronic sector. 

4.4.1 Effective hadronic weak mixing angles sin^ d^^ 

For the four lightest quarks (u, d, c, s) the SM calculation is as advanced as in the leptonic 
sector (see eq. (62)) [16]. For the bottom sector the calculation is more involved, as an 
additional mass scale enters the calculation due to top-quark dependent two-loop vertex 
graphs. The resulting additional leading top mass dependent two-loop terms can be accounted 
for via [73, 75] 



where is the form factor of the down-quark. The contribution parametrises the difference 
between and k^. It is given by [73,75] 




(65) 




167r2 



(67) 
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The explicit expression for ' can be found in Ref. [73,75]. 

In our MSSM calculation we include the available higher-order SM and MSSM contribu- 
tions in complete analogy to Sect. 4.3.1, see in particular eq. (63)^. In the hadronic process 
Z qq, SUSY QCD contributions enter at the one- loop level already (see Fig. 18). As 
explained above, we include these corrections into the form factors. Numerically the SUSY 
QCD corrections only play a subleading role. The Zbb vertex graphs with virtual Higgs ex- 
change contain couplings that are enhanced by tan (5. We resum the leading contributions as 
described in Sect. 2.3. We have furthermore included the two-loop SM contributions to Kb of 
eq. (67) (the one-loop terms contained in eq. (67) have been subtracted from the full MSSM 
one-loop result to avoid double-counting). 



4.4.2 Hadronic decay widths 

As a first step to calculate the hadronic partial widths we derive pq from the effective one-loop 
couplings. The incorporation of SUSY QCD corrections, the resummation of tan /3-enhanced 
contributions and the inclusion of further higher-order corrections proceeds in the same way 
as described above. Leading non-universal corrections to the ZU) vertex are obtained with [73, 
75] 

P5 = Pd(l+T-5)^ (68) 

where again parametrises the difference between down- and bottom-quark couplings as 
in eq. (67). 

The radiation factors Ry are more involved than in the leptonic case, since they incor- 
porate both final state QED and QCD interactions, as well as the bottom-quark mass in the 
process Z bb. For light quarks with ~ the form factors are of comparatively simple 
form [54] 

v,A ^ ^ y ^ J y ^ J -^?4 ^2 

(69) 

The explicit expressions for Ryj^, including also the case rUq ^ 0, were implemented as given 
in Refs. [4,54]. 

A further correction to the Z boson width available in the literature arises from non- 
factorisable two- loop contributions. The mixed non-factorisable QCD and electroweak cor- 
rections for u,d,c,s-quarks are taken from Ref. [22]. For the bottom-quark we use results 
from Ref. [23]. The non-factorisable corrections can simply be added to the respective par- 
tial hadronic width. 



^For the u,c,d- and s-quarks only parametrisation formulas for sin^ arc available in the literature. We 
therefore extract Ak^ ' from these formulas by applying the same strategy as in Ref. [30]. 
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4.5 Decay width for Z X^X? 

As discussed above, as final ingredient for tlie computation of the total Z boson width in the 
MSSM we evaluate the partial width for the process Z This is done in analogy to 

the leptonic case described above. For we use 

which is in accordance with the expression in eq. (65) for = 0. 



5 Numerical analysis 

We now present our numerical results for the W boson mass, M^^, and the most relevant 
Z boson observables: the effective leptonic weak mixing angle, sin^6'eff, the total Z boson 
width, Vzi the ratios for the leptonic and 6-quark width of the Z boson, Ri and i?;,, and the 
hadronic peak cross section, o"°ad- explicitly discuss the effective hadronic weak 

mixing angles defined in Sect. 4.4.1, which nevertheless enter the hadronic decay widths via 
the quantities k/, see eq. (36) ^ . We start our numerical discussion with a detailed investi- 
gation of the predictions for M^^^ sin^ ^eff, T^, i?;, 7?^, and a'^^ with respect to the different 
SUSY masses and complex phases. Then we discuss the results for these observables in the 
MSSM for a choice of sample scenarios and discuss their decoupling behaviour with respect 
to the SM limit. Specific scenarios such as the CPX scenario [35] and "Split SUSY" [33] are 
investigated. Finally a scan over all relevant SUSY parameters is performed. Also the SUSY 
contributions to the invisible Z boson width are analysed. 

The numerical analysis of our analytical results for the Z boson observables, which were 
calculated as described above, and M\\r, sec Rcf. [30], is performed with the help of a newly 
developed Fortran program called SUSY-POPE (SUSY Precision Observables Precisely Eval- 
uated), which will be made publicly available [31]. Though built up from scratch, for the 
calculation of the MSSM particle spectrum our code partially relies on routines which are part 
of the FormCalc [61] package. The Higgs sector parameters are obtained from the program 
FeynHiggs [40,43,44]. 

If not stated otherwise, in the numerical analysis below for simplicity we choose all soft 
SUSY-breaking parameters in the diagonal entries of the sfermion mass matrices, eq. (2), to 
be the same, 

Mf = Mp^Mp,^ ... . (71) 
In the chargino/ neutralino sector the GUT relation 

Ml = \'\m2 (72) 

(for real values of Mi and M2) is often used to reduce the number of free MSSM parameters. 
We have kept Mi as a free parameter in our analytical calculations, but will use the GUT 
relation to specify Mi for our numerical analysis if not stated otherwise. 

^All Z observables evaluated in this paper, including also forward-backward and left-right asymmetries, 
have already been included in a recent analysis in the constrained MSSM (CMSSM) [12]. 
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We have fixed the SM input parameters as 



= 1.16637 X 10"^ 

1/137.03599911, 
mt = 170.9 GeV [79] , 

For the bottom-quark mass, m^, we list the pole mass given above as a scheme-independent 
reference point. Following Ref. [4], in the numerical calculations we use the corresponding 
running bottom quark mass mi,{Mz) ~ 2.8 GeV. 

The results for physical observables are affected only by certain combinations of the com- 
plex phases of the parameters the trilinear couplings A^, A^, . . . , and the gaugino mass 
parameters Mi, M2, M3 [80,81]. It is possible, for instance, to rotate the phase away. 
Experimental constraints on the (combinations of) complex phases arise in particular from 
their contributions to electric dipole moments of heavy quarks [82], of the electron and the 
neutron (see Refs. [83,84] and references therein), and of deuteron [85]. While SM contribu- 
tions enter only at the three-loop level, due to its complex phases the MSSM can contribute 
already at one-loop order. Large phases in the first two generations of (s)fermions can only be 
accommodated if these generations are assumed to be very heavy [86] or large cancellations 
occur [87], see however the discussion in Ref. [88]. Accordingly (using the convention that 
4>M2 = 0- done in this paper), in particular the phase 0^ is tightly constrained [89], while 
the bounds on the phases of the third generation trilinear couplings are much weaker. 

5.1 MSSM parameter dependence 

We start by comparing our full MSSM result for the EWPO {M^, sin^Oes, Tz, Ri, Rb and 
"^had) ^i^^ current experimental results, which are hsted in Tab. 1. In the following we 
present results where the most important SM and SUSY parameters are varied in order to 
identify the observables with the highest sensitivity to SUSY loop effects. 

5.1.1 Dependence on the sfermion mass scale 

In Fig. 3 we show the prediction for the EWPO for real parameters as a function of and 
indicate how this prediction changes if the top-quark mass is varied within its experimental 
la interval, mt = (170.9 ± 1.8) GeV [79]. The other parameters are A,6,t = 2Mj, fi — 
Ma — mg — M2 — 300 GeV and tan /3 = 10. The result is compared with the current 
experimental values, see Tab. 1. It can be seen in Fig. 3 that only Mw, sin^ ^eflf and Tz 
exhibit a pronounced sensitivity to Mj. While Mw shows a mild preference for light M^, 
sin^ ^^cflf is in better agreement with the experimental value for large Mj. In interpreting 
the latter result it should be noted that the world-average on sin^ ^^cff involves independent 
measurements that differ from each other by more than three standard deviations [1]. An 
experimental resolution of this issue will most likely require an ILC with a GigaZ option. 
The prediction for Tz lies within its observed la error for most of the parameter space. The 
dependence of Ri, Rt,, and a^^^ on Mj is nearly flat, where the first two observables are 
within the la band, and the latter is slightly below. The impact of varying mt within its 



Mz = 91.1875 GeV, a,(Mz) = 0.118, 

a42i = 0.02758 [76, 77], Aaiep = 0.031498 [78], 
m6 = 4.7GeV, = 1.777 GeV, 

mr — ms — . . . — 0. 
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observable 


central exp. value 


^ — ^today 


^LHC 
(7 


^ILC 
(7 


Mw [GeV] 


80.398 


0.025 


0.015 


0.007 


sm e/eflf 


0.23153 


0.00016 


0.00020-0.00014 


0.000013 


Tz [GeV] 


2.4952 


0.0023 




0.001 


p 


on 7fi7 


u.uzo 




n ni 

U.Ul 




0.21629 


0.00066 




0.00014 


"^had 


41.540 


0.037 




0.025 



Table 1: Summary of the clcctrowcak precision obscrvablcs that will be analysed in the 
following, including their current experimental central values and experimental errors, a = 
(jtoday Also showu are the anticipated experimental accuracies at the LHC, a^^^ and 

at the ILC (including the GigaZ option), cx^^^. Each number represents the combined results 
of all detectors and channels at a given collider, taking into account correlated systematic 
uncertainties, see Refs. [6,90-92] for details. A recent review can be found in Ref. [93]. Non- 
existing analyses are referred to as " — " . 

experimental error is non-negligible only for Mw, sin^ ^eff, and Fz- It results in the following 
shifts^ 

5mr = 1.8GeV [79] ^sin" = 5.4 x 10-^ (75) 

Similarly we have analysed the impact of a shift in Aa^^i) 

^^para,AaSd ^ 5 3 ^eV, (77) 

5( Aajfij) = 3.5 X 10-^ [1, 76, 77] 5 sin^ ^P^^^-^^Sa ^ 12 x 10-^ (78) 

^j^para,Aa,^d = 0.32 MeV. (79) 

It should be noted that the parametric uncertainty in sin^ ^eff induced by (5(Aq;^^'|j) is of similar 
size as the current experimental uncertainty. A significant improvement of the uncertainty 
in Aq;[^^^^ is clearly very desirable in order to be able to fully exploit future progress in the 
experimental measurements of the EWPO as well as in their theoretical predictions, see also 
Sect. 6. 

*The parametric uncertainties were evaluated for the SPSla' benchmark point, see Sect. 5.3, to allow for 
a comparison with the theoretical errors estimated in Sect. 6. However, the uncertainties are only weakly 
dependent on the mass scale of the SUSY particles. 
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Figure 3: Prediction for Mw, sin^6'efr, Tz, Ri, Rb and cr^a^j as function of the common 
sfermion mass scale Mj for mt = (170.9 ± 1.8) GeV. The SUSY parameters are tan/5 = 10, 
Ar = At = Ab = 2M?, Ma = l^ = mg = M2 = 300 GeV. 
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5.1.2 Dependence on fj, and M2 

The variation of the EWPO with the parameters from the chargino/neutrahno sector is 
investigated in Fig. 4. We show the prediction for the six EWPO as a function of M2 for 
Ijl = 250,500, 1000 GeV (and Mi is chosen according to the GUT relation). The other 
parameters arc set to Mf = 300 GeV, At^b^^ = '^Mj, Ma = 1000 GeV, tan/5 = 10 and 
rrig = 600 GeV. As expected, also in this case the observables with the largest sensitivity to 
the variation of the SUSY parameters are Mw, sin^ 9es and Tz- The impact of varying M2 
and fjL on the other three observables is neghgible. The variation with fj, results in shifts in 
MiY, sin^^eff and at the la level. A sizable variation with M2 can only be observed for 
M2 ^ 200 GeV for this set of parameters. While M]y (sin^ ^eflf) shows a monotonous decrease 
(increase) by ~ lo" with M2, Tz exhibits a strong increase up to M2 < 150 GeV and then 
slowly decreases for further increasing M2. Good agreement between the Tz prediction and 
the experimental value is found for small M2, see also Sect. 5.6. 

5.1.3 Dependence on complex phases 

For the analysis of the dependence of the complex phases we focus on the three EWPO 
that show the strongest variation with the SUSY parameters, M^, sin^ ^eff and Tz, see the 
two previous subsections. Since the dependence on the sfermion mass parameters is much 
stronger than on the chargino/higgsino parameters we only investigate the dependence on the 
phases of At and A^. As for Ar/Myy [30] we find that the effective one-loop couplings g^y^} (1) 

depend only on the absolute values \Xh\ of the off-diagonal entries in the t and b mass 
matrices, where = A^ — n/ tan/3, Xh = A^ — fj, tan/?. Thus, the phases of /x, A^ and 
enter only in the combinations (0^^ + <f>n), giving rise to modifications of the squark masses 
and mixing angles. It furthermore follows that the impact of (pAt (0a J on the sfermion masses 
(see eq. (6)) is stronger for low (high) tan/5. 

In Fig. 5 we show the three EWPO as a function of 4>At (with 0^ = (pAf, = 0, left plots) 
and (f)A^ (with 0^ = (f)At = 0, right plots) for different values of tan P (varied from tan (5 — h 
to tan/3 = 45). The other parameters are set to — Mh± — M2 — rrig — 500 GeV, 
l^t.fe.rl = Ia*| = 1000 GeV, = = 03 = 0. As expected, the dependence of My/, 
sin^ ^cfj and Tz on (pAt is most pronounced for small tan/5 (see left panel of Fig. 5). The 
variation of (pAt in this case gives rise to a shift in the three precision observables by 1-2(7. 
The effect becomes smaller for increasing tan /5, up to tan [3 = 15. On the other hand, for high 
tan /3 the lighter b mass becomes rather small for the parameters chosen in Fig. 5, reaching 
values as low as about 100 GeV for tan/5 = 45. This leads to a sizable shift of ~ 1-2(7 
in the EWPO already for vanishing phases. The slight rise in the dependence on (pAt for 
tan /5 > 25 is due to the overall enlarged SUSY contributions which occur for large tan /5 and 
the resulting low sbottom masses. 

The dependence of the precision observables on (plots on the right-hand side of Fig. 5) 
is rather small, except for the highest tan/5 value shown in Fig. 5, tan/5 = 45. This is again 
related to the sizable correction induced by the lighter b mass. In this scenario the variation 
of 4>At, yields a shift in Mw, sin^ 6'eff and Tz at the la level. 
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Figure 4: Prediction for Mw, sin^6'efr, ^z, Ru Rb and a^^^^ as function of M2 for /i = 
250, 500, 1000 GeV. The remaining SUSY parameters are tan j3 = 10, Mj = 300 GeV, 
Ar = At = Ab = 2Mp m-g = 600 GeV, Ma = 1000 GeV. 
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5.2 Impact of Z 



As next step in our numerical analysis we investigate the impact of the decay Z — » XiXi 
on the invisible Z boson width, Finv- A sizable ZxiXi coupling in combination with a light 
neutralino, m^o ^ Mz/2, results in an additional contribution to the invisible width of the 
Z boson, in addition to the SM decays into neutrinos. Since the experimental result for 
the invisible width of the Z boson is somewhat below the SM prediction [1], additional 
contributions from new physics are tightly constrained. It is of interest in how far the 
possibility of a light neutralino is affected by the precision measurement of the invisible 
width of the Z boson. 

The mass of the lightest neutralino is determined by Mi, M2 and fi, see eq. (10) (here 
we assume all parameters to be real). If Mi or M2 or /z is much smaller than the other two, 
the lightest neutralino is mostly a bino, zino, or higgsino, respectively. Its mass is to a large 
extent determined by this smallest mass value. If the condition 

tan (3 M2 



Ml = Mf ^ = 2M|s 



2 „2 



1 + tan2 /3 ) (/iM2 - 2M| tan/3/(l + tan2 /5 )c2 



0) 



was exactly fulfilled, the lightest neutralino, Xi? would be massless and almost entirely bino- 
hke. 
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Figure 6: The invisible Z boson width Finv in- and excluding the process Z XiXi is 
shown as a function of Mi (left) and rriyXi (right). The SUSY parameters are chosen to be: 
M^- = 250 GeV, = At = A^, = = m-g = 500 GeV, Ma = 500 GeV, M2 = 200 GeV. A 

light, almost entirely bino-like, Xi is obtained by varying Mi around the value m[^^ defined 
in eq. (80) within the range of -100 to +100 GeV. 

In Fig. 6 we analyse a scenario with a light neutralino (see also the discussion in Ref. [94]). 
We do this by varying Mi around the value m{°-* defined in eq. (80), i.e. 

Ml = Mf ^ + AM, (81) 
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Figure 7: The invisible Z boson width Finv in- and excluding the process Z XiXi is 
shown as a function of Mi (left) and m^o (right). The SUSY parameters are chosen as 
follows: Mj = Ar = At = Ab = Ma = 'm-g = 600 GeV, ^ = 125 GeV, Ma = 200 GeV, 
Ml = -100... 100 GeV. 

where we let AM run from —100 to +100 GeV. The other parameters in Fig. 6 are Mj = 
250 GeV, At,b,T = /i = = Ma = 500 GeV, Ms = 200 GeV, tan/3 = 10. In the left (right) 
plot of Fig. 6 the results for Fjnv are shown as a function of Mi (m^o). The curves in the 
plot on the right hand side of Fig. 6 consist of two branches for a given value of m^o. This 
is a consequence of the fact that different values of Mi in general lead to different values 
of Finv although the value of m^o may be the same. Comparing with the plot on the left 
one can see that the upper branch corresponds to Mi ^ 0, the lower branch to Mi ^ 0. 
The coupling of a nearly pure bino to Z bosons is very weak, resulting in a visible, but 
negligible effect on Finv- The entire range of Finv shown in Fig. 6 is outside its la error. As 
mentioned above, this is a consequence of the well-known fact that the measured value of the 
invisible Z width is below the SM prediction [1]. As the predictions including and excluding 
the process Z XiXi differ at most by ~ 0.02 MeV, compared to an experimental range 
of = 499.0 ± 1.5 MeV, the inclusion of the decay into neutralinos Xi clearly does not 
give rise to additional constraints on Mi or the mass of the lightest neutralino m^o in this 
scenario. 

A different scenario with /i ~ Ml, but Ml < iMa (in coarse agreement with the GUT 
relation) is analysed in Fig. 7. The parameters are set to M^ = At^b,T = Ma = rrig = 600 GeV, 
H = 125 GeV, Ma = 200 GeV, tan/? = 10, and Mi is varied from -100 GeV to +100 GeV. 
In this scenario the lightest neutralino can have a sizable higgsino component. In the left 
(right) plot of Fig. 7 the results for Finv are shown as a function of Mi {m^o). The lower 
branches in Fig. 7 (plot on the right) again correspond to Mi ^ 0, the upper branches to 
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Ml > 0. For values m^o < Mz/2 the F^o contribution to Finv becomes sizable. In this 
case the deviation between the MSSM prediction for the invisible Z width and the measured 
value, which is slightly above la if the decay into neutralinos is not open, raises above the 
2(7 level for |Mi| < 35 GeV or m^o < 35 GeV. 

As a consequence, if the SUSY parameters are such that Xi has a sizable coupling to 
the Z boson (e.g. M2 > Mi ^ /x), the theoretical prediction for Finv can easily exceed the 
experimental central value by several standard deviations (provided that m^o < Mz/2). In 
such a scenario the contribution to the invisible width arising from F^o yields interesting 
constraints on the parameters /i, Mi, and M2. 

5.3 The SPS benchmark scenarios 

In this section wc analyse the six EWPO Mw, sin^ ^eff, F^, Ri, R^, and a^^^ in the SPS la', 
SPS lb and SPS 5 benchmark scenarios [32,95]. Our analysis extends the results of Ref. [30] 
to the Z pole observables. In order to analyse the dependence of the EWPO on the scale of 
supcrsymmctry we scale for each SPS point all SUSY parameters carrying mass dimension 
by a common factor, i.e. Ma — (scalefactor) x M^^, Mp p, — (scalefactor) x M|^|,, A,6,t = 
(scalefactor) x A^^^, /i = (scalefactor) x /i^^^, Mi_2,3 = (scalefactor) x Mff|. In Fig. 8 we 
show Mw, sin^^eff, ^z, Ri, Rb, and a^^^ as a function of the lighter t mass, nif^, in the 
three SPS scenarios. As before, only Mw, sin^ 6'cff and Tz show a sizable variation with m^^ 
(i.e., the scalefactor). For all the EWPO the (loop-induced) variation between the three SPS 
scenarios is relatively small. Mw shows best agreement with the experimental results for low 
rrif^. sin^ ^eff lies in the ±1(7 range for mf^ ^ 400 GeV. F^ shows a variation with mf^ of 
~ 1 a around the experimental value. 

Fig. 9 shows for the example of the SPS la' scenario that, as expected, the MSSM predic- 
tions for the EWPO approach the corresponding predictions in the SM (for M^ = M^^^^) 
for large values of the SUSY mass scale. The predictions within the MSSM and the SM 
(for M^ — Mjf^^^) are shown as a function of m^^ (as before, all parameters carrying 
mass dimension are scaled by a common factor). The variation of the SM prediction with 
the SUSY scale is induced by the corresponding change in the light CP-even Higgs-boson 
mass. The difference between the MSSM and the SM predictions becomes negligible for 
mf^ ^ 2000 GeV for Mw and sin^ ^eflf- For the other EWPO the decoupling of the supersym- 
metric contributions occurs for even lower values of . For the hadronic peak cross section, 
(T^g^^, the change with the SUSY scale in the MSSM and the SM prediction has opposite signs 
(on the other hand, for the individual factors entering o-h^d' ^1- (^2), the dependence of 
the MSSM and SM predictions on the SUSY scale goes into the same direction). The fact 
that we recover the most up-to-date SM prediction for the EWPO in the decoupling limit is 
a consequence of the procedure described in Sect. 4.2.1 (see eq. (59)). 
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Figure 9: Predictions for Mw, sin^6'efr, Tz, Ri, Rb and a^^^ within the SPSIa' scenario in 
comparison with the SM result calculated for = Mh- The observables are shown as a 
function of m^^, the mass of the lighter of the two scalar top-quarks. The SPS parameters of 
mass dimension are varied with the scale of supersymmetry as described in the text. 
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5.4 The CPX benchmark scenario 



In this section we analyse the prediction of the EWPO in the CPX scenario. This scenario 
is defined as [35] 

= 500 GcV, \At\ = 1000 GcV, = A = 
M2 = 200 GeV, \mg\ = 1000 GeV, ^ = 2000 GeV (82) 
^At,i.,r = 09 = 7r/2 , 

with the aim to indicate the possible size of CP-violating effects. 

The LEP Higgs searches have been interpreted within the CPX scenario [96]. A parameter 
region for intermediate tan/5 and and light M^^ could not be excluded at the 95% C.L., so 
that no lower hmit on the mass of the hghtest Higgs boson could be set. The reason for 
the appearance of these "CPX holes" is a strong suppression of the coupling of the lightest 
Higgs boson to gauge bosons, while the second lightest Higgs boson (also having a somewhat 
reduced coupling to gauge bosons), which may also be within the kinematic reach of LEP, 
decays predominantly into the lightest Higgs boson. The latter leads to a rather difficult 
topology of the final state. 

As a consequence of its strongly suppressed coupling to gauge bosons, the impact of 
the contributions of a light Higgs boson associated with the "CPX holes" on the EWPO is 
expected to be rather small. In Fig. 10 we investigate whether the current measurements of 
the EWPO, see Tab. 1, allow to put constraints on the parameter space of the "CPX holes". 
The consistent inclusion of the loop-corrected Higgs boson masses and coupfings is crucial 
for this analysis. 

In Fig. 10 we show the predictions of the six EWPO as a band for the range tan ^ — 3- 
13, where the blue curves (boxes) correspond to tan/? = 3, the purple curves (crosses) 
to tan P = 13. The band is to the left cut off for the lowest physically allowed values 
of Mh± ■ (For illustrative purposes we also show the values of the three neutral Higgs boson 
masses in the lowest plot of Fig. 10 as obtained with FeynHiggs [40,43,44].) The regions 
unexcludcd by the LEP Higgs searches (starting at tan/5 ^ 3.6) are indicated by a light blue 
shading for the two observables Mw and sin^ ^^cfr (for simplicity, we have omitted them for 
the other observables). As explained above, our one- loop result contains the full complex 
phase dependence, whereas at the two-loop level we use the approximation formula given in 
eq. (60). While this method works for tan/3 = 13, for tan/3 = 3 the value ior (f) — n (see 
eq. (60)) cannot be evaluated since (at least) one of the squark mass squares turns negative. 
Therefore for tan/5 = 3 we use 0(0) = 01^(0) + (©^"^^(O) - O^^{0)), O = Mw, sin^^cfr, 
Tz, Ri, Rb, o"had- This results in a higher theoretical uncertainty of the prediction for the 
tan/5 = 3 border shown in Fig. 10, which should be kept in mind for the interpretation of 
the figure. 

Fig. 10 shows that the predictions for the EWPO in this parameter region of the CPX 
scenario are in general in good agreement with the experimental results. Only for low tan /5 
(i.e., values close to the boundary of the region indicated by the contour with tan/5 = 3) 
and small Mh± sizable deviations can be observed, most notably for sin^ ^^gfr- On the other 
hand, the parameter regions corresponding to the "CPX holes" show only small deviations 
from the EWPO. Most of the corresponding regions are even within the 1 a intervals of the 
experimental values of Mw and sin^ ^eff- Thus, more precise EWPO measurements combined 
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case are shown. 
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with improved theoretical predictions (and correspondingly smaller intrinsic uncertainties in 
the EWPO calculations) would be needed to reach the sensitivity for probing the "CPX hole" 
regions via their effects on EWPO^. 

5.5 Scenario where no SUSY particles are observed at the LHC 

It is interesting to investigate whether the high accuracy achievable at the GigaZ option of 
the ILC would provide sensitivity to indirect effects of SUSY particles even in a scenario 
where the (strongly interacting) superpartners are so heavy that they escape detection at the 
LHC. 

We consider in this context a scenario with very heavy squarks and a very heavy gluino. 
It is based on SPS la', but the squark and gluino mass parameters are fixed to 6 times their 
SPS la' values. The other masses are scaled with a common scale factor as described in 
Sect. 5.3, except M4 which we keep fixed at its SPS la' value. In this scenario the strongly 
interacting particles are too heavy to be detected at the LHC, while, depending on the 
scale-factor, some colour-neutral particles may be in the ILC reach (the reach for the direct 
production of colour-neutral particles at the LHC will be rather limited [97]). In Fig. 11 we 
show the prediction for sin^^cfr in this SPS la' inspired scenario as a function of the lighter 
chargino mass, m-±. The prediction includes the parametric uncertainty, o-p^ra-iLC^ induced 
by the ILC measurement of rrit, Snit = 100 MeV [98,99], and the numerically more relevant 
prospective future uncertainty on Aaj^'^j, 6{Aal^^^) = 5x 10~^ [100]. The MSSM prediction for 
sin^ ^efif is compared with the experimental resolution with GigaZ precision, a^^'^ — 0.000013, 
using for simplicity the current experimental central value. The SM prediction (with M^^ — 
^MSSM^ is also shown, applying again the parametric uncertainty o-p^ra-iLC^ 

Despite the fact that no coloured SUSY particles would be observed at the LHC in this 
scenario, the ILC with its high-precision measurement of sin^ 9es in the GigaZ mode could 
resolve indirect effects of SUSY up to m^± < 500 GeV. This means that the high-precision 
measurements at the ILC with GigaZ option could be sensitive to indirect effects of SUSY 
even in a scenario where SUSY particles have neither been directly detected at the LHC nor 
the first phase of the ILC with a centre of mass energy of up to 500 GeV. 

5.6 Heavy scalar masses 

The scenario discussed in the previous section is characterised by rather heavy scalar quarks 
(and a heavy gluino). Other scenarios with heavy scalar masses that found attention in 
recent years are the so-called "split SUSY" scenario [33] and the "focus point region" [34] of 
the CMSSM. For completeness, we also briefly discuss the sensitivity of the EWPO to these 
scenarios. 

In the split SUSY scenario the fermionic masses (i.e. the chargino, neutralino, and gluino 
masses) are relatively light and retain their GUT-induced hierarchy, resulting in a heavier 

^In case that a sizable deviation between the predictions in the CPX scenario and the measiirements of 
the EWPO would occur, it would be important to check whether changes in the CPX parameters (e.g. shifts 
in the mass parameters of the first two families) could bring the EWPO prediction into agreement with the 
experimental results, while not aff'ecting the LEP Higgs analyses (which are mostly affected by the parameters 
of the third family). 
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Figure 11: Theoretical prediction for sin^ ^eff in the SM and the MSSM (including prospective 
parametric theoretical uncertainties) compared to the experimental precision at the ILC with 
GigaZ option. An SPS la' inspired scenario is used, where the squark and gluino mass 
parameters are fixed to 6 times their SPS la' values. The other mass parameters are varied 
with a common scalefactor, see Sect. 5.3. 



gluino and lighter charginos and neutralinos. The scalar mass parameters, on the other hand, 
are set to very large values, larger than 10^-10^ GeV, i.e., they decouple from the predictions 
of the EWPO. Consequently, only a small deviation in the EWPO prediction from the SM 
limit is to be expected in this scenario. We focus here on the two most sensitive EWPO, Mw 
and sin^ ^eff (the results given for Mw are an update of Ref. [30]). 

In Fig. 12 we show in a split SUSY scenario the SUSY contribution, i.e., the difference 
between the SUSY result and the corresponding SM result, to Myy (upper plot) and to sin^ 9es 
(lower plot). This is done by choosing a large value for Mj and subtracting the SM result 
with Mfj^ = M^s^^ from the result obtained in the MSSM. We have chosen Mj = 5 TeV 
and Ma = 2.5 TeV. Choosing even higher values for the scalar mass parameters would only 
lead to negligible shifts as compared to the results shown Fig. 12. The results are displayed in 
Fig. 12 in the /X-M2 plane for tan (3 = 10. The gluino mass has been fixed to rrig = 500 GeV 
(the results are insensitive to this choice). The results are given in terms of contour lines 
representing the different collider precisions, see Tab. 1. The region excluded by chargino 
searches at LEP [58] is inside the dashed white lines. As can be seen in the figure, deviations 
from the SM prediction of more than 30 MeV in Mw and 3 x 10"'* in sin^ ^cfr occur only 
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Figure 12: 6Mw = (M^^sm _ ^^^^^^ pj^^) dsin^Oes = (sin^ ^,^|ssm _ sin^ el^) 

(lower plot) in a split SUSY scenario. The SUSY parameters are = 5000 GeV, Ma = 
2500 GeV, tan /3 = 10, = At = Ab = 2Mj, m-g = 500 GeV. Contour lines are drawn 
for 6Mw = {30 MeV, a^^^ = 15 MeV, a^^^ = 7 MeV, 3 MeV, 2 MeV} and ^sin^^es = 
{30 X 10-^ 15 X 10~^ 3(T^L^ = 3.9 x 10^^ 2(tIlc = 2.6 x 10-^ a^^^ = 1.3 x 10^^}. 



36 



in those experimentally excluded areas. The parameter regions leading to a shift as large 
as the GigaZ precisions in Mw of 1(7^^° = 7 MeV and in sin^^eff of la^^^ = 0.000013 are 
significantly larger. In particular, a 1 a effect in sin^^gfr can occur for IM2I ^ 400 GeV 
(for negative (//M2) only). (Similar results were obtained in Ref. [101], see also Ref. [102].) 

Another scenario with heavy scalar masses is the focus point region [34] in the CMSSM 
(the CMSSM is characterised by a common scalar mass parameter, mo, a common fermionic 
mass parameter, mi/2, and a common trilinear coupling, Aq, at the GUT scale, supplemented 
by the low-scale parameter tan/3 and the sign of the parameter //). In the focus point 
region mi/2 is relatively small, while itlq is rather large, itlq ~ few TeV, and also tan/? is 
relatively large, tan/? ^ 40. We have investigated (using the program ISAJET 7.71 [103]) 
the predictions for M^, sin^ ^^cfr and Vz in the focus point region. The results for sin^ ^eff 
and Vz follow the pattern observed in Ref. [30], where My^ had been analysed. The main 
contributions to the shifts in the EWPO arise from the chargino and neutrahno sector (see 
Sect. 5.1.2) and only very small effects arise from the scalar fermion contributions. For 
tan/3 = 50, > 0, mi/2 = 250 GeV, mo = 1500 GeV, Aq = —250 GeV, corresponding to 
a point with the currently lowest value of mi/2 in the focus point region for which the dark 
matter density is allowed by WMAP and other cosmological data (see, for example, Ref. [9] 
for a more detailed discussion), we find an effect of the SUSY contribution of 6 x 10~^ in 
sin^ ^eff and 1 MeV in Fz- Thus, even for the low-mi/2 region of the focus point scenario 
the GigaZ precision for sin^ ^eflf will be needed to gain sensitivity to indirect effects in the 
EWPO. 

5.7 Higgs sector at higher orders 

We next investigate the impact of higher-order contributions associated with the Higgs sector 
of the MSSM. As mentioned above, since the Higgs sector enters the EWPO only via loop 
corrections, in order to evaluate predictions for the EWPO at one-loop order it would formally 
be sufficient to treat the Higgs sector in leading order, i.e., at tree level. However, at the 
tree level the predicted value for the mass of the light CP-even Higgs boson is so low that 
it is below the exclusion limit from the Higgs searches at LEP [96, 104]. As a consequence, 
treating the MSSM Higgs sector at tree level in the predictions for the EWPO would lead to 
artificially large contributions to the EWPO from the light MSSM Higgs boson. Therefore 
a consistent incorporation of higher-order contributions in the Higgs sector, as described in 
Sect. 2.3, is crucial in order to be able to use realistic mass values for the light MSSM Higgs 
boson and to take into account potentially large higher-order effects. 

To analyse the numerical effects of higher-order contributions in the Higgs sector we 
study Mw, sin^6'eff, Tz, and Rf, for the MSSM parameter set tan/3 = 50, M2 = 300 GeV, 
m-g = 600 GeV, Mf = 300 GeV, // = 300 GeV, \At\ = \Ah\ = \A^\ = 2M;, = tt/2, 
fpAr = 4>At, = 4>tJ. = 0Mi = 4>g = 0. Mh± is varied from 85 GeV to 1000 GeV. A rather large 
value for tan (3 was chosen to further analyse a possible numerical impact of the resummation 
of leading tan/3 enhanced Higgs (s)bottom couplings in the process Z bb, entering via 
Aft, see eq. (16). Our results are shown in Fig. 13, where "Higgs Born" labels the numerical 
results where only Born-level Higgs masses and couplings were used. "Higgs Full" are the 
results which take into account the implementation described in Sect. 2.3, i.e., they account 
for loop-corrected Higgs masses and mixing angles, mixing of CP eigenstates in the presence 
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of complex MSSM parameters, and tan (5 enhanced Higgs (s)bottom couplings, included via 

In the upper row of Fig. 13 we show the results for Mw and sin^ 6'efr- The impact of 
the higher-order corrections in the Higgs sector, corresponding to the difference between the 
two results shown in Fig. 13, amounts to a shift in and sin^ 6'eff of about la. Thus, 
this contribution, which is formally of two-loop order, has a sizable numerical impact and 
should be taken into account in order to arrive at a precise prediction for Mw and sin^ ^eff- 
The lower row of Fig. 13 displays the results for Vz and i??,, where we also show the impact 
of the A5 corrections alone by comparing with the result where = (labeled as "Higgs 
Full, Afe = 0"). The numerical impact of the higher-order corrected Higgs boson sector on 
the observables Vz and i?^ relative to their current experimental errors is less pronounced 
as compared to Mw and sin^^eff- Setting A;, = yields only a relatively small shift from 
the full result. It should be kept in mind that Af,, see eq. (16), approximates the leading 
contribution from the sbottom sector only for large SUSY mass scales. If the relevant particle 
masses are simultaneously small, < mj, the theoretical uncertainties in the predictions 
of the EWPO can be slightly larger. 

5.8 MSSM parameter scans 

As a final step of our numerical analysis we investigate the behaviour of the two EWPO that 
are most sensitive to higher-order effects in the MSSM, Mw and sin^ ^efr, by scanning over 
a broad range of the SUSY parameter space. The following SUSY parameters are varied 
independently of each other in a random parameter scan within the given range: 

sleptons : Mp p, = 100 . . . 2000 GeV, 
hght squarks : Mp = 100 . . . 2000 GeV, 

' up / down 

t/h doublet : Mp p, = 100 . . . 2000 GeV, 

' up/down 

/l^^^ = -2000 . . . 2000 GeV, 
gauginos : Mi,2 = 100 . . . 2000 GeV, 
= 195 . . . 1500 GeV, 
II = -2000 . . . 2000 GeV, 
Higgs : Ma = 90 . . . 1000 GeV, 

tan /3 = 1.1... 60. (83) 

Those parameters which may in general be complex are taken to be real, as the dominant 
effects of the complex parameters only enter via the shifts they induce in the sparticle masses 
and mixings (see the discussion in Subsect. 5.1.3). Effects of this kind are thus covered 
by scanning over the parameters and ranges given cq. (83). Performing the scans, only the 
constraints on the MSSM parameter space from the LEP Higgs searches [96,104] and the lower 
bounds on the SUSY particle masses from direct searches as given in Rcf. [58] were taken 
into account. Apart from these constraints no other restrictions on the MSSM parameter 
space were made. 

In Fig. 14 we compare the SM and the MSSM predictions for sin^ ^eff as a function of 
rrit as obtained from the scatter data. The predictions within the two models give rise to 
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Figure 13: Predictions for Mw, sin^ ^efr, ^z, and Rb as a function of Mh± . The MSSM param- 
eters are tan/5 = 50, Ma = 300 GeV, m-g = 600 GeV, = 300 GeV, /i = 300 GeV, \At\ = 
\Ab\ = \Ar\ = 2Mj, = vr/2, = 0a, = 0^ = 0^ = 0, Mh± = 85 GeV . . . 1000 GeV. 
"Higgs Born" labels the results calculated for Born-level Higgs sector parameters, while 
"Higgs Full" refers to the results calculated for the loop-corrected Higgs sector parameters 
as described in Sect. 2.3. For Tz and Rb (second row) also the result without resummation 
of tan P enhanced contributions to the bottom Yukawa coupling is shown, labeled as "Higgs 
Full, Ab = 0". 



two bands in the m^-sin^ 6'eff plane with only a relatively small overlap region (indicated 
by a dark-shaded (blue) area). The allowed parameter region in the SM (the medium- 
shaded (red) and dark-shaded (blue) bands) arises from varying the only free parameter of 
the model, the mass of the SM Higgs boson, from Mf^ = 114 GeV, the LEP exclusion 
bound [104] (lower edge of the dark-shaded (blue) area), to 400 GeV (upper edge of the 
medium-shaded (red) area). The very light-shaded (green), the light shaded (green) and the 
dark-shaded (blue) areas indicate allowed regions for the unconstrained MSSM. In the very 
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Figure 14: MSSM parameter scan for sin^ O^s as a function of nit over the ranges given in 
eq. (83). Todays 68% CL. ellipses as well as future precisions, drawn around todays central 
value, are indicated in the plot. 
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Figure 15: MSSM random parameter scan for Mw as a function of rrit over the ranges given in 
eq. (83). Todays 68% CL. ellipses as well as future precisions, drawn around todays central 
value, are indicated in the plot. 
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Figure 16: MSSM random parameter scan over the ranges given in eq. (83). The top- 
quark mass is treated as an additional free parameter and is varied in the range rrit = 
165 . . . 175 GeV. Shown is the combination of Mw and sin^ 9es- Todays 68% C.L. elhpses as 
well as future precisions, drawn around todays central value, are indicated in the plot. 

light-shaded region at least one of the ratios mi^/mi^ or Tni^/mi^ exceeds 2.5,^ while the 
decoupling limit with SUSY masses of 0{2 TeV) yields the upper edge of the dark-shaded 
(blue) area. Thus, the overlap region between the predictions of the two models corresponds 
in the SM to the region where the Higgs boson is light, i.e., in the MSSM allowed region 
{Mh ^ 130 GeV [43,44]). In the MSSM it corresponds to the case where all superpartners 
are heavy, i.e., the decoupling region of the MSSM. The 68% C.L. experimental results 
for rrit and sin^ 6'efr are indicated in the plot. As can be seen from Fig. 14, the current 
experimental 68% C.L. region for rrit and sin^ 6'eff is in good agreement with both models 
and does not indicate a preference for one of the two models. The prospective accuracies 
for the Tevatron/LHC (5 sin^ ^2™*™"/^"^ = 0.00016, ^^Tevatron/LHC ^ ^ ^^^^ j^C 

with GigaZ option (5 sin^ 6'^^^^'^^^^^ = 0.000013, drn^'^^'^^^^'^ = 0.1 GeV) are also shown in 
the plot (using the current central values), indicating the strong potential for a significant 
improvement of the sensitivity of the electroweak precision tests [91]. 

In Fig. 15 we compare the SM and the MSSM predictions for Mw as a function of rrit as 
obtained from the scatter data''. The ranges of the varied parameters and the band structure 
of the SM and MSSM predictions are analogous to the ones in Fig. 14. The experimental 
value for Mw includes the latest CDF measurement [106], resulting in the world average 



^ We work in the convention that m; < m; . 

Jl — }2 

The plot shown here is an update of Refs. [6,30, 105]. 



41 



of Mw = (80.398 ± 0.025) GeV [2, 106]. The 68% C.L. region for and Mw exhibits a 
preference for the MSSM over the SM. The prospective accuracies for the Tevatron/LHC 
((5mJ^'™"'°''/^"^ = 15 McV) and the ILC with GigaZ option {5M^^'^'^'^^'^ = 7 MeV) are also 
shown in the plot (using the current central value). 

Finally in Fig. 16 we show the combination of Mw and sin^ with the top-quark mass 
varied in the range of 165 GeV to 175 GeV. The ranges of the other varied parameters and 
the colour coding are the same as in Figs. 14, 15. The current 68% C.L. experimental results 
for Mw and sin^ ^eflf are indicated in the plot. The region of the SM prediction inside todays 
68% C.L. ellipse corresponds to relatively large rrit values, outside the current experimental 
range of rrit = (170.9 ± 1.8) GeV [79]. Thus, the combination of Mw and sin^ 6'eff exhibits 
a slight preference for the MSSM over the SM. Again also shown are the anticipated future 
improvements in the measurements of Mw and sin^ ^eff- 

6 Remaining higher-order uncertainties for sin^ ^eff 

Following the discussion in Ref. [30] we now estimate the missing higher order uncertainties 
in the prediction of sin^^^cff- As seen in Sect. 5, besides the W boson mass, sin^ ^^cfi is the 
observable with the most pronounced dependence on higher-order SUSY contributions. Thus, 
it is important to reduce its intrinsic theoretical uncertainty from unknown higher-order (SM 
and SUSY) correctins sufficiently below its experimental error and its parametric theoretical 
uncertainties (induced by the experimental errors of the input parameters). The remaining 
SM theory uncertainty was estimated to be [13, 16] 

5sin2^f^ = 4.7x 10-^ (84) 

This corresponds to the theory uncertainty of our prediction for sin^ 6cs in the decoupling 
limit of the MSSM, since we have incorporated all known SM higher order contributions into 
our result (see Sect. 4.3). 

As detailed in Refs. [28,30], additional theoretical uncertainties arise from higher-order 
corrections involving supersymmetric particles in the loops. Depending on the overall SUSY 
scale these uncertainties were estimated to be [28] 

5sin2 ^eff = 4.7 X 10-^ for Mf < 500 GeV, 

5 sin^ ^eff = 1.5 X 10-^ for Mf = 500 GeV, (85) 
^sin^ ^eff = 1.3 X 10-^ for = 1000 GeV. 

Adding SM and SUSY uncertainties in quadrature one finds 5sin^^eff — (4.9 — 6.6) x 10~^, 
depending on the SUSY mass scale [28]. 

Additional theory uncertainties arise for complex parameters, as we only include the 
full phase dependence at the one-loop level. However, MSSM two-loop terms which are 
only known for real MSSM parameters are incorporated into the full prediction for sin^ ^eff 
via the interpolation relation eq. (60). Following the prescription in Ref. [30], where the 
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corresponding error in the M\y evaluation had been obtained, we estimate here the maximal 
error over the £ [0, tt] interval to be^ 

5sin2 e^s = 2.6 X 10"^ for < 500 GcV, 

6 sin^ ^cff = 1.7 X 10"^ for Mj = 500 GeV, (86) 
dsin^ ^eff = 0.6 X 10"^ for Mj = 1000 GeV. 

Similar values are found in an independent approach: the approximation formula eq. (60) can 
be applied to the one-loop case, see Sect. 4.2.3. The difference between the approximation 
and the full phase dependence at one-loop order is scaled to the two-loop level by applying 
a conservative factor of 0.2 (obtained from the analyses of the SUSY contributions to the 
p parameter [26,28]). The result for the estimate of the uncertainty induced by the approxi- 
mation formula for the phase dependence at the two-loop level is similar (slightly below) to 
the numbers in eq. (86). 

The full theoretical uncertainty from unknown higher-order corrections in the MSSM with 
complex parameters can now be obtained by adding in quadrature the SM uncertainties from 
eq. (84), the theory uncertainties from eq. (85) and the additional SUSY uncertainties from 
eq. (86). This yields 6sm^ O^s = (4.9 - 7.1) x 10~^ depending on the SUSY mass scale. 

The other source of theoretical uncertainties, besides the one from unknown higher-order 
corrections, is the parametric uncertainty induced by the experimental errors of the input 
parameters. The corresponding uncertainties for 5mt = 1.8 GeV and ^i^ct^^^) are given in 
eqs. (74) - (79). The uncertainty in rrit will decrease during the next years as a consequence 
of a further improvement of the accuracies at the Tevatron and the LHC. Ultimately it will be 
reduced by more than an order of magnitude at the ILC [98,99] , see also the phenomenological 
analysis in Rcf. [107]. For A«|j^j one can hope for an improvement down to 5 x 10~^ [100], 
reducing the parametric uncertainty to the 1.8 x 10^^ level (for a discussion of the parametric 
uncertainties induced by the other SM input parameters see, for example, Ref. [6]). In order 
to reduce the theoretical uncertainties from unknown higher-order corrections to the level of 
1.8 x 10~^, further results on SM-type corrections beyond two-loop order and higher-order 
corrections involving supersymmetric particles will be necessary. 



7 Conclusions 

We have presented the currently most accurate evaluation of Z pole observables in the MSSM. 
These comprise the effective weak mixing angle, sin^6'efr, Z decay widths to SM fermions, 
r{Z — > //), the invisible and total width, Finv and F^, forward-backward and left-right 
asymmetries, Aps and Alr, and the total hadronic cross section, o"had- The calculation 
includes the complete one-loop results, for the first time taking into account the full complex 
phase dependence. It furthermore includes all existing higher-order corrections available 
in the MSSM. We have also incorporated all available SM contributions that go beyond 

^As representative SUSY scenarios we have chosen SPS la', SPS lb, and SPS 5, each for Mf = 1000 GeV, 
500 GeV, and for Mf < 500 GeV. The lowest values considered for Mf are roughly 300, 300, 400 GeV for 
SPS la', SPS lb, SPS 5, respectively. For lower values the parameter points arc excluded by Higgs mass 
constraints. The light stop mass for the SPS 5 point lies considerably below 400 GeV. 
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the existing MSSM corrections. These corrections can be numerically important since the 
evaluations of higher-order corrections in the SM are more advanced than in the MSSM. As 
a result, in the decoupling limit (i.e., if all SUSY mass parameters are large) our calculation 
reproduces the currently most up-to-date SM results for the Z pole observables. Concerning 
the Z boson decay widths we obtained for the first time a full one-loop result for V{Z X?X?) 
that contributes to the invisible width of the Z boson if this decay is kinematically possible. 
A pubhc computer code based on the results for electroweak precision observables (the Z pole 
observables and the W boson mass [30]) is in preparation [31]. 

We presented a detailed numerical analysis of the impact of the various MSSM sectors 
on the prediction of the EWPO. We find that the EWPO M^y, sin^ Q^<s and V z are sensitive 
to variations of the Ijh sector parameters and to a lesser extent of the chargino/neutralino 
parameters. The impact on the other Z pole observables was found to be much smaller 
and well within one experimental standard deviation. The evaluation of and sin^ ^cfr 
incorporating higher-order corrections to Higgs boson masses and couplings differs by about 
1 a from the predictions where the tree-level approximation for the Higgs sector is used. 

Particular emphasis in the analysis was put on the impact of the complex phases at the 
one-loop level. The largest phase dependence arises in general in the scalar top sector. Large 
effects are also possible in the scalar bottom sector if tan /3 is large. Shifts induced by a phase 
variation exceeding one experimental standard deviations can occur for M14/, sin^ ^cfr and F^, 
depending on the other MSSM parameters. It should be noted, however, that the dependence 
on the phases appearing in the lib sector at the one- loop level drop out of the explicit one- loop 
terms for the effective electroweak couplings S'^y^j (j^) and thus only enter via their impact on 
the squark masses and mixing angles. Consequently, a measurement of the CT'-conserving 
EWPO alone cannot serve to reveal the presence of the complex phases. However, precise 
measurements of the EWPO in combination with other experimental information will be very 
valuable for constraining the SUSY parameters including the complex phases. 

By comparing our MSSM result with the corresponding SM result we have illustrated the 
sensitivity of the EWPO to the virtual effects of SUSY particles. We pointed out that the 
anticipated experimental accuracy at future colliders, in particular at the ILC with GigaZ 
option, might resolve the virtual effects of SUSY particles, even in scenarios where the SUSY 
particles are so heavy that they escape direct detection at the LHC and the first phase of the 
ILC. 

The EWPO have been analysed in the framework of several MSSM scenarios. As typical 
examples we investigated three SPS-like scenarios. Further studied scenarios were the "focus 
point" region, "split SUSY" and the CPX scenario. In the "focus point" and "split SUSY" 
scenarios, both containing heavy scalar fermions, only very small deviations from the corre- 
sponding SM result could be observed. Within the CPX scenario we found that the current 
experimental precision of the EWPO is not yet sufficient to probe the parameter space at 
low Higgs boson masses that could not be covered by the LEP Higgs searches. 

Finally, we have analysed the theoretical uncertainty in the sin^ Q^<g^ prediction arising from 
the incomplete inclusion of complex phases beyond the one-loop level. Our estimate yields 
that this (additional) uncertainty can amount up to 2.6 x 10~^, depending on the SUSY 
parameters. The combination of this new uncertainty with the one estimated previously 
for the MSSM with real parameters yields an estimate of the intrinsic uncertainty of up to 
7 X 10~^ for small SUSY masses. 
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The results derived for the purpose of this paper have already found application in various 
precision analysis projects [9,10,12], which are aimed at providing indirect constraints on 
the scale of supersymmetry based on current constraints from cosmology, B physics, and 
electroweak precision data. 
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A Appendix 

Here we show generic Feynman-diagrams contributing to the Z pole observables in the MSSM. 
A.l One-loop vertex graphs for the process Z — ^ // 




Figure 17: Generic vertex contributions to Z ^ ff. 
A. 2 SUSY QCD one-loop graphs for the process Z ^ qq 




Figure 18: Generic SUSY QCD vertex contributions to Z — > qq. 
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A. 3 One- loop graphs for Z — x^Xi 
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Figure 19: Generic vertex contributions to Z — > x?x?. 



A. 4 SUSY two- loop graphs contributing to Ap 
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Figure 20: Sample diagrams for the SUSY 0{aas) contributions to Ap: [g] squark loop with 
gluon exchange, {g) (s)quark loop with gluino exchange. 
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Figure 21: Sample diagrams for the three classes of MSSM 0{ol\)^ 0{o?^)^ 0{atab) contribu- 
tions to Ap: (q) quark loop with Higgs exchange, (q) squark loop with Higgs exchange, (H) 
quark/squark loop with Higgsino exchange, (f denotes Higgs and Goldstone boson exchange. 
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